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We study in this work a scenario that the universe undergoes a two step phase transition with the first step 
happened to the dark matter sector and the second step being the transition between the dark matter and the 
electroweak vacuums, where the barrier between the two vacuums, that is necessary for a strongly first order 
electroweak phase transition (EWPT) as required by the electroweak baryogenesis mechanism, arises at the 
tree-level. We illustrate this idea by working with the standard model (SM) augmented by a scalar singlet 
dark matter and an extra scalar singlet which mixes with the SM Higgs boson. We study the conditions for 
such pattern of phase transition to occur and especially for the strongly first order EWPT to take place, as 
well as its compatibility with the basic requirements of a successful dark matter, such as observed relic density 
and constraints of direct detections. We further explore the discovery possibility of this pattern EWPT by 
searching for the gravitational waves generated during this process in spaced based interferometer, by showing 
a representative benchmark point of the parameter space that the generated gravitational waves fall within the 


sensitivity of eLISA, DECIGO and BBO. 


PACS numbers: 11.30.Er, 11.30.Fs, 11.30.Hv, 12.60.Fr, 31.30.jp 


I. INTRODUCTION 


The Standard Model (SM) of particle physics so far pro- 
vides an excellent description of a wide variety of experi- 
mental observations, nevertheless, it must be extended be- 
cause it fails at least in explaining two cosmological puz- 
zles: it does not provide a dark matter (DM) candidate and 
it is not possible to generate the observed baryon asymme- 
try of the universe (BAU). The later one requires a strongly 
first order electroweak phase transition (EWPT) to provide a 
non-equilibrium environment if the BAU is generated via the 
electroweak baryogenesis mechanism (Hah. Since the 125 
GeV Higgs boson is too heavy to give rise to a first order 
EWPT is a new ingredients beyond the SM Higgs inter- 
actions may be needed. On the other hand, how DM interacts 
with the SM particles remains unknown. The Higgs portal 
turns to be attractive since it may kill two birds (DM and the 
strongly first order EWPT (7) with one stone. 


The characteristics of a scalar DM triggered strongly first 
order EWPT are twofold: (i) It is a two-step phase transi- 
tion (sif14} and the barrier between electroweak symmetric 
and broken phases arises at the tree-level, which may avoid the 
problem of gauge-dependence [15]; (ii) The parameter space 
of the model is strongly constrained by the exclusion limits of 
direct detections and can not explain the observed relic den- 
sity. It was pointed out in Ref. that, for the SM plus a real 
scalar DM S scenario, S could only constitute up to 3% of 
the total DM relic density, while still be relevant for the direct 
detection because of its sizable coupling with the SM Higgs 
as required by the two-step EWPT. 
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In this paper, we revisit this type of EWPT by extending 
the SM with a scalar singlet DM S and another scalar singlet 
® which can mix with the SM-like Higgs boson. Due to the 
mixing of ® with the SM Higgs, there are two separate con- 
tributions to DM nucleon scattering, which may cancel with 
each other and leads to a negligible cross section. It results in 
a mechanism to suppress the direct detection signal and evade 
the currently most stringent experimental constraints from the 
direct detection experiments ira (iat. Furthermore, the quar- 
tic interaction 6?,S? may contribute to the mass of S, which 
can lead to a relatively small quartic coupling of S with the 
SM-like Higgs and result in a sizable relic density. As a result 
the tension between the DM relic density and direct detec- 
tion which occurs for the SM plus singlet case can be 
highly loosed in this scenario. 


For the thermal history of the universe, the presence of the 
two additional scalars allows a two step EWPT. This hap- 
pens as follows: As the universe cools down, S gets vac- 
uum expectation value (VEV) first and the universe transits 
to this phase; As the temperature gets lower, a second mini- 
mum occurs in the SM-like Higgs and ® subspace. The uni- 
verse then tunnels to this electroweak minimum resulting in 
a strongly first order EWPT with gravitational waves (GW) 
generated. GW signals coming from this two-step EWPT are 
calculated using one benchmark point of the parameter space 
and our results show that the corresponding GW signals are 
testable in the Evolved Laser Interferometer Space Antenna 
(eLISA), DECi-hertz Interferometer Gravitational wave Ob- 
servatory (DECIGO), Ultimate-DECIGO and Big Bang Ob- 
server (BBO). 


The remaining of this paper is organized as follows. Con- 
ventions of the model are defined in Sec. We then study 
the DM phenomenology in Sec. and the capability of a 
strongly first order EWPT in Sec. The GW calculations 
and discovery prospects are explored in Sec.|V| after which 
we present a brief summary in Sec. [VI] 
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Il. THE MODEL 


We present the model for two-step phase transition and 
scalar DM in this section. The model extends the SM with 
two real scalar singlets: S and ®, where S' is the DM candi- 
date, while ® gets non-zero VEV at the zero temperature and 
thus mixes with the SM Higgs (68). The Higgs potential can 
be written as 


1 1 1 1 
Vv, ye 2 p2 z pt _ p2 Q2 = 4 
0 aa + 7% 5HsS + 7rs8 
—wWHIH + d(HtH)? +),S°H'H + °H H 
+3576", (1) 


where H is the SM Higgs doublet. The effective potential, 
critical for the EWPT, can be written as: Veg = Vo + Vow + 
Vr, where Vo is the tree-level potential, Vcw is known as 
the Coleman-Weinberg term Gil, Va includes finite tempera- 
ture contributions from loops |22] and bosonic ring PABA. 
For one-step phase transition, where Higgs potential contains 
no tree-level cubic term, the barrier between the electroweak 
symmetric and broken phases usually arises from loop correc- 
tions and according to the Nielsen identity B5), the resulting 
effective potential is gauge dependent. For the two step phase 
transition, the barrier arises at the tree-level. So we include 
in the effective potential the standard one-loop T # 0 cor- 
rections, but only retain terms proportional to O(T") for the 
consideration of gauge invariance. Thermal masses of scalars 
take the following forms: 


o 3g? + 9°? À ht Ay + A2 2 
m = { ie aai w fi: 
As à AB] me 
i242 4 ale 
{+ 3 =| (3) 
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e A E ee 4 
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We require S to have zero VEV at T = 0 for the stabil- 
ity of the DM and parametrize the other two fields by H = 
(0, (vew +h)/V2) and ġ = va +¢ with here vew ~% 246 GeV 
and h; ~ 1. Minimization conditions around (vew, va, 0) in 
the (h, œ, s) space allow us to trade the two VEVs vew and va 
for u? and u4 by 
H? = Anvew +A2VB, Hg = A2vEW + ABS - (5) 
We also replace 42, by the physical DM mass m2, 
L = Atiy + 2A30% — m2. (6) 


With these substitutions, the mass matrix for (h,@) is then 
given by 


M2 =2( Yew 
VEwUo A2 


which can be diagonalized by a 2 x 2 orthogonal matrix 
parametrized by a rotation angle 0. The mass eigenstate (h, 8) 
can then be written as 


h= coh+ so, 


VEwUo A2 ) (7) 


vere 


ĝ = —soh + cod. (8) 


2 


where cg = cos@ and sg = sin@. We define h as the SM- 

like Higgs boson and therefore when 0 — 0, h is purely the 

SM Higgs. With this definition, three of the parameters in 

the potential A, A2, Ag can be replaced by the physical masses 
2 


2 : 
m$, ms and 6: 


22 2 42 
Me Co + 4M 86 


cs 9 
wey ; (9) 

m? — m2 
àz = sin(29) > 2, (10) 

AVEWVe 

m? 83 + mcs 

lo = L1, 11 
d TA (11) 


In summary, we have a total of 7 free parameters: 


ms, As; 0, M1; A3, (12) 


Uð, m 3 
where A, and Àg are relevant for the DM phenomenology. 

Before proceeding to study EWPT and DM phenomenol- 
ogy, we discuss constraints on the mixing angle 0 from 
Higgs measurements and electroweak precision measure- 
ments. Couplings of the SM-like Higgs to all SM particles 
are rescaled by the factor cg due to the mixing. As a result, 
signal rates, the ratio of Higgs measurements relative to the 
SM-like Higgs measurements, equal to c?. Performing a uni- 
versal Higgs fit to the data given by the ALTAS and CMS 
collaborations, one has || < 0.526 at the 95% CL [26]. It 
was shown in Ref. that the constraint of electroweak pre- 
cision is closely related to the mass of db, which is weaker than 
that of Higgs measurements for a light b and turns to be much 
strong as mg gets heavier. For constraints of vacuum stabil- 
ity and perturbativity, we refer the reader to Refs. for 
detail. 


HOI. DARK MATTER 


The fact that about 26.8% of the universe is made of DM, 
whose relic abundance is about Qh? = 0.1189 (30), has been 
well established, while the nature of the DM is still unknown. 
The Higgs potential in Eq. (I), has Z2 x Zə discrete symme- 
try for S and ® respectively. At the zero temperature, ® gets 
non-zero VEV and thus its corresponding Z2 is broken. On 
the other hand, we require that the global minimum at T = 0 
forbids a nonzero VEV in the s direction. As a result s can 
serve as a DM candidate, which interacts with the SM parti- 
cles via the Higgs portal interactions. The couplings between 
the DM and the physical scalars are 


s? h : 2A3ug80 + A1vewco , (13) 

so : 2A3U6CO — AL VEWS¢ , (14) 
A 1 

sh? : zu + A353 ; (15) 
> 1 

ep o: 5186 + A3ce i (16) 


s*h $ : cgse(2A3 — A1) . (17) 
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FIG. 1: Results from a scan over the parameter space in plane (ms, A1) for the left panel and (0, A1) for the right panel. Every point in these 
plots generates the desired EWPT pattern and leads to vanishing DM nucleon scattering cross section at tree level. The green points further 
gives DM relic density within the interval Qh? € (0.03, 0.12). The red points, in addition to generating relic density falling into this interval, 


also satisfy the strong first order EWPT condition v, (Toc)/To = 1. 


The relic density of the DM is governed by the Boltzmann 
equation: 


n+ 3Hn = — (ov) (n? — nko) , (18) 


where H is the Hubble constant, (ov) is the thermal average 
of the reduced annihilation cross section, For our model, the 
DM can annihilate into f f, ZZ/WW, hh, ġġ and hę final 
states via s-channel Higgs mediations, where di-scalar pro- 
cesses also receive contribution from four point interactions in 
Eq. . Compared with the most simple Higgs portal, here 
there are more annihilation channels for a heavy DM. Cou- 
plings in Eqs. and are both relevant to the DM direct 
detections, while the contribution of b mediated processes is 
suppressed by the factor s3 and thus is negligible for a small 
0. For numerical simulations we implement the model in Lan- 
HEP and use MicrOMEGAs to calculate the 
relic density. 

For direct detections, the spin independent scattering cross 
section off nucleon is given by 

2 


2,2 F 
pom Coa; Soa 2 T 
on =- aad F 9 + 9 5 fr, | a9) 
TURwits | mM a q=u,d,s 


where aj,, a3 are effective couplings in Eqs. and 
respectively and ÍT, are nucleon form factors for light 
quarks BS). The inclusion of ® that mixes with the SM Higgs 
can lead to a complete cancellation in contributions to on at 
tree level. This occurs when the effective coupling of DM 
with quarks vanishes, that is, when the expression in |---| 
vanishes. Imposing this condition would allow us to eliminate 
one more parameter A3, 


r UpwA1 (m? tan + m3 cot 0) ne 
3 =— Qve (m? = m=) , ( ) 


which will be adopted in the following numerical analysis for 
simplicity. 

For a detailed survey of the parameter space, we perform 
a parameter scan and show the results in Fig. [I] in the plane 
(ms, Àı) (left panel) and in the plane (8, A1) (right panel). 
In these plots, the green points give the relic density within 
the range Q-h? € (0.03, 0.12) allowed by the current results 
from Planck Bd, in addition to being able to induce a two- 
step EWPT pattern as discussed later. From the plot in the left 
panel, we can see the coupling À; that can give a relic density 
within this range is of magnitude (0.005, 0.02) for a relatively 
heavy DM, while it can be larger ~ 0.06 for a relatively light 
DM. For the right plot, the mixing angle is mostly negative 
since most of the scanned points are for m g > Ms, and there- 
fore a negative 0 is needed to give a positive \2 according to 
Eq.[10] The red points in both plots also give a strongly first 
order EWPT condition up,(Tc)/Te 2 1. 


IV. ELECTROWEAK PHASE TRANSITION 


The barrier between the electroweak symmetric and broken 
phases emerges at the tree-level in our two-step phase tran- 
sition scenario. As a result, the effective potential, which in- 
cludes the standard one-loop T ¥ 0 corrections but only retain 
the leading terms( thermal mass terms) in the high T' expan- 
sion so as to avoid problems relating to the gauge dependence, 
can be written as 

l, 2 2 1 4 2 
Veg = 3 (He Ig) 3 (Hs II; )s 


1 1 1 1 
= 2 TI, )h? + rau + seh? + zA h? 


1 i 
+A387¢? + pre + 7A, (21) 
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FIG. 2: This figure gives an illustrative picture of the two step EWPT 
with the first step in the s direction and the subsequent one from s 
direction to the (h, d) direction. 


in terms of background fields h, ¢ and s. 


Given Eq. (21), one can trace the evolution of the universe 
phase as temperature drops and our desired pattern of EWPT 
is illustrated in Fig.|2} At sufficiently high temperature, the 
universe sits at the global minimum (0, 0,0) where the elec- 
troweak symmetry is restored. As temperature drops to T,, a 
minimum develops in the s direction and due to an absence 
of a barrier with the one at the origin, the universe transits to 
this minimum through a second order phase transition. As T 
continues decreasing to Th, a second minimum develops in 
the (h, d) direction but its free energy is initially higher than 
the one in the s direction. As T further decreases, the free en- 
ergy of the minimum in the (h, d) direction drops faster than 
the one in the s direction, and at the critical temperature Tc, 
these two minima become degenerate. Slightly below Tc, the 
universe makes a second transition to the global minimum in 
the (h, ġ) direction. Due to the existence of the barrier be- 
tween these two minima, the transition may occur as a first 
order EWPT by a tunneling process and proceeds through nu- 
cleations of electroweak bubbles which expand, col- 
lide and coalesce, leaving eventually the universe in the elec- 
troweak broken phase. To avoid washing out of the generated 
baryons by the sphaleron process inside the electroweak bub- 
ble, the sphaleron process needs to be sufficiently suppressed 
inside the bubble and this translates into the following gener- 
ally adopted criterion ABBA: 


vn(To) 
To 


>1. (22) 


Note that due to the reflection symmetry of the effective po- 
tential under h + —h, 6 + —ġ or s — —s, there are identical 
phase structures and we focus on the region with h, ¢, s > 0. 


The VEV of the SM Higgs at the finite temperature can be 


written as 


Uh (T) = (23) 


The critical temperature Tç can be calculated analytically 
from the following equation: 


Alua — Tg)” — 2A2(u? — Mn) (ua — He) + àa (u? — IMa)? 
Me- X 


where ju”, vz, and u% can be written in terms of physical pa- 
rameters. 

As a concrete example, we show in the left panel of Fig. B] 
the evolution of V at these two minima as the tempera- 
ture drops from right to the left, for the parameter choice 
ve = 65GeV, mz = 82GeV, ms = 71GeV, As = 0.015, 
0 = 0.12, Ay = 0.046 and A3 = 0.57. Here the horizontal 
green line on the far right denotes the symmetric phase at high 
temperature when the universe sits at the origin, the magenta 
dashed line represents the minimum in the s direction and the 
blue line is the minimum in the (h, œ) direction which even- 
tually evolves to the electroweak minimum at the zero tem- 
perature labeled by a red dot. The s direction phase appears 
at Ts = 150GeV continuously away from the origin while 
the (h, ġ) direction phase starts from (159GeV, 9.15GeV, 0) 
at Tn = 108GeV above which it is a saddle point. The T at 
which the blue and magenta lines intersect is the critical tem- 
perature To = 78.8GeV and is labeled by a black dot. At Tc, 
the two minima are (205GeV, 45GeV, 0) and (0, 0, 352GeV). 
For more details on the evolution of the minimum in the (R, ¢) 
direction, we plot in the right panel the track of this mini- 
mum in blue dotted line for the whole thermal history with 
here the red arrow denoting the direction as temperature drops. 
The contours in this plot gives an measure of the value V at 
T = 0 and the electroweak vacuum labelled by the red point 
sits at the deepest location. For this parameter choice, the 
strongly first order EWPT criterion Eq.(22) is achieved since 
un(To)/To = 2.6. 

A more comprehensive survey of the model is given by a 
scan over the parameter space as shown in Fig. [L| where all 
plotted points give above described pattern of EWPT by im- 
posing various conditions during the scan. These conditions 
include: (i) there are two minima in the field space with one 
in the s direction and the other in (h, ¢) direction; (ii) the 
electroweak minimum needs to be lower than the one in the 
s direction at T = 0 should the s direction minimum persist 
at T = 0; (iii) the minimum in the s direction occurs earlier 
than the one in the (h, ) direction, that is, Ts > Th. Af- 
ter imposing these conditions, the points are further filtered to 
give a relic density in the range (0.03, 0.12) with the remain- 
ing plotted using green color. At this step, these points can 
give a two-step EWPT and a sizable amount of relic density 
but the second step EWPT is not necessarily strongly first or- 
der. As such, we further calculate the critical temperature Tce 
and vp (Tc) corresponding to the black intersection point of 
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FIG. 3: The left panel shows the evolution of V at the two minima in s(blue) and in (h, ¢)(magenta dashed) directions as T drops from right 
to left. The right panel shows the tracks of the minimum (blue dotted line) in the (h, œ) direction for the thermal history with the contours 
denoting the values of V at T = 0. In these plots, the magenta, blue, black and red dots represent the temperatures at Ts = 150GeV, 
Ty = 108GeV, To = 78.8GeV and T = 0 respectively. The parameters are chosen as: ve = 65GeV, m = 82GeV, ms = 71GeV, 
As = 0.015, 6 = 0.12, Ar = 0.046 and A3 = 0.57 which gives v, (Tc)/Te = 2.6. 


the blue and magenta lines for each above green point. The 
points which satisfy the condition in Eq. are shown with 
red color in Fig.[I]. So we can see there are sufficient param- 
eter space in this model where the desired strongly first order 
EWPT pattern can be realized. We note in passing that the 
second order phase transition for the first step could also be 
strongly first order, should a more complete effective poten- 
tial be adopted. 


V. GRAVITATIONAL WAVES 


The calculation of the GWs generated during the second 
step EWPT needs a numerical analysis of the tunneling pro- 
cess at finite temperature and in particular involves solving 
the critical bubble profiles. As mentioned in previous sec- 
tion, after the universe cools down to a temperature below 
Tc, the second transition from (0, 0, vs) to the true vac- 
uum (vh, Ug, 0) takes place by the nucleation of true vacuum 
bubbles. This tunneling rate per unit time per unit volume 
reads 


T ~ A(T)e—53/7, (24) 


where S3 is the Euclidean action of the critical bubble which 
minimizes the action 


> 2 
S3(6,7) = 4r | Pdr (42) +V(G.7)] . 05 


and the prefactor A(T) is roughly of O(T*), whose pre- 
cise evaluation needs integrating out fluctuations around the 


bounce solution il. From this rate formula, the bubble nu- 
cleation temperature T» is defined as the probability for a sin- 
gle bubble to be nucleated within one horizon volume being 
O(1), i.e. 


A oo 4 
J Tvd = | a (=) e Ps = O(1),26) 
0 Tn 


where Vy (t) is the horizon volume, Mp; is the Planck mass 
and Ç ~ 3 x 1072.. This equation implies numerically 
S3(Tn)/Tn œ~ 140 and serves as our definition of T),. 
From the bounce solutions, two other important parameters œ 
and £, that are directly relevant for the calculation of GWs, 
are defined by: 


== ’ B = A, Th, a ’ 
dT le 


"ae 
Prad 


(27) 


where Pyac is the vacuum energy density released in the 
phase transition, H,, is the Hubble parameter at Tn, p74 = 
gxm*T+/30 with g, the relativistic degrees of freedom in the 
plasma at T,,. A small 3/H,, will trigger strong phase transi- 
tion and consequently a significant stochastic background of 
GWs. 

The observable characterizing the GW background is the 
energy spectrum h?Qgw(f) and is given by 


B h? dpgw 
pe dlog f’ 


where pew is the energy density of GWs with f its frequency 
and pe is the critical energy density today. The sources of 
stochastic GW signals arising from first order EWPT can be 
classified into three categories: (1) collisions of bubble walls 


h?Qew(f) 


(28) 
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FIG. 4: The left panel shows changes of S(T’) /T in the neighborhood of T, and the right panel shows the GWs generated during the first order 
EWPT as a function of frequency from three sources: sound waves(blue dotdashed line), collision(red dotted line), turbulence(green dashed 
line) and total contribution(cyan solid line). The color shaded regions fall within the experimental sensitivities of eLISA(four configurations 
with notation NiAjMkLI), ALIA(gray), BBO(green), DECIGO(yellow) and Ultimate-DECIGO(purple). 


and shocks in the plasma; (2) sound waves in the plasma af- 
ter the bubble collision; (3) Magnetohydrodynamic (MHD) 
turbulence in the plasma [43]. The total energy spectrum is 
given approximately by the sum of these three sources, 


h?Qaw od h?Qeon + ome T h? Desai j (29) 


The GW contribution from bubble collisions can be calculated 
using the envelope approximation [4446], which numerically 
results in the following GW contribution to the spectrum (471: 


HO? / ka \? /100\3 
h? Quo = 1.67 x 1075 | — 
nls) = na (Fe) G) (FP) 


0.1103, 3.8(f/feou)?® 
i E + =) f + el s] i 


where v,, is bubble wall velocity, x characterizes the frac- 
tion of latent heat deposited in a thin shell and feon is 
the peak frequency produced by the bubble collisions. At 
the time of the phase transition fy = 0.628/(1.8 — 
0.luw + v2,), which is red-shifted to give the peak fre- 
quency today: feon = fy x a(Tn)/ao = 16.5 x 
1076 Hae" fe yy 100 GeV) (ga (100)! Both v, and 
8] 


coll 
k are functions of a, which read as | 


_ 1/V3 + Va? + 2a/3 


l+a 


mw OTa + 0.278 
~ 1+0.7a 


The sound wave contribution of the GW intensity is numeri- 
cally fitted by 


; 7 Le (Hn \ / sva \? (100\ 
Path) = 265 10° (He) (E2 (10) 


FA? 7 7/2 
me (4) Tage p ea 


31) 


w 


where ky denotes the faction of latent heat transformed into 
bulk motion of the fluid, fsw is the peak frequency that can 
be given from the rescaling of its value at the phase transition, 
ie. fow = fay X a(Ta)/a0 = (2/V3)(8/vw) X a(Tn)/ao. 
We refer the reader to Refs. for the value of k, in the 
small and large vw limit. The MHD turbulence contribution 
to the GW spectrum can be written as 


P T af A g 
h2Qturn(f) = 3.35 x 10 ( : ) (=) r 
— f 
(1+ f/f) H(A + 8rf/hn) ? 


where ktu © 0.lku (50), hn is the Hubble parameter today, 
and ftu ~ (3.5/2)(8/vw) x a(Tn)/ao. 

We use the package CosmoTransitions to solve for pro- 
files of the critical bubble and the nucleation temperature for 
the benchmark point used earlier. This parameter choice gives 
Q.h? = 0.04 constituting about 34% of the total DM relic 
density. For this benchmark parameter point, we show in the 
left panel of Fig. [4]the quantity $(T)/T as a function of T 
and we found T, = 41.2GeV. The resulting GW signals are 
shown in the right panel of Fig.[4]where the three contributions 
are shown with details in the caption. As can be seen from 
this figure, the blue sound wave contribution dominates and 
is almost indistinguishable from the sum of the three sources 
which is denoted by the cyan line. For a comparison with ex- 
periment, we show firstly four sensitive regions correspond- 
ing to four configurations of the eLISA detector. These four 
regions are labelled in the format NiAjMKkLI and are plotted 
as the red shaded regions at the top following conventions of 
Ref. (431 52). We can see for this parameter set, the GWs can 
be detected by the configuration N2A5M5L6 while unreach- 
able by the others. We further add the sensitivitive regions of 


X Uw (33) 
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several other proposed decihertz GW experiments: Advanced 
Laser Interferometer Antenna (ALIA) (53 (69), BBO, DE- 
CIGO and Ultimate-DECIGO . The data are taken from 
Ref. and are plotted as gray, green, yellow and purple 
regions respectively . From the plot, we can see the GW 
signals from the benchmark point fall within the detectable 
ranges of BBO, DECIGO and Ultimate-DECIGO. It should 
be noted that there might be other points in the model param- 
eter space where the generated GWs can have much larger 
energy density spectrum and are therefore within reach of the 
other three eLISA configurations as well as ALIA, and give 
a sufficiently large DM relic density in the meanwhile. This 
however needs a dedicated survey of the parameter space in- 
cluding all the considerations and we leave it to future works. 
For GWs from one-step EWPT in SM extended with scalar 
singlet(s) il, we refer the reader to Refs. as well 
as references cited in these papers for detail. Alternatively, 
first order phase transition may lead to originations of primor- 
dial black holes. They can be captured by neutron stars or 
astrophysical black holes [sd], resulting in GWs that may be 
detected by Advanced LIGO or Advanced Virgo. 
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